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' '  Time-domaln  dipole  antenna  surray  studies  conducted  at  the  Electromagnetic  Effects  i 
Laboratory^  of  the  Harry  Diamond  Laboratories  for  application  to  large  scale  nuclear  electro¬ 
magnetic  pulse  (EMP)  simulation  are  described.  Analytically  and  experimentally  derived  EMP 
environments  for  a  single-dipole  antenna  are  compared  with  corresponding  data  for  multiple 
element  broadside  and  collinear  dipole  antenna  arrays.  The  multiple-element  concept  is 
studied  as  a  technique  for  elimiimting  the  need  for  extremely  high  voltages  in  the  single  ele¬ 
ment-type  nuclear  EMP  simulators  by  substituting  the  superimposed  electromagnetic  fields 
radiated  by  several  lower  voltage  sources  to  obtain  the  equivalent  EMP  environment.  Experi¬ 
mentally  measured  input  impedances  are  compared  with  the  results  from  classical  antenna 
theory  to  determine  the  effects  of  mutual  coupling  between  elements  in  the  array  configura¬ 
tions.  The  aperture  size  of  the  array  is  shown  to  be  an  important  factor  in  determining  the 
pulse  shape  characteristics  of  the  EMP  environment,  in  the  near  and  far  fields.  Results  of 
the  studies  define  specific  advantages  and  limitations  in  the  dipole  antenna  array  technic^ies 
for  nuclear  EMP  simulation. 
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ABSTRACT 

Time  domain  dipole  antenna  array  studies  conducted  at  the  Electromag¬ 
netic  Effects  Laboratory^  of  the  Harry  Diamond  Laboratories  for  application 
to  large  scale  nuclear  electromagnetic  pulse  (EMP)  simulation  are  described. 
Analytically  'Jid  experimentally  derived  EMP  environments  for  a  single 
dipole  antenna  are  compared  with  corresponding  data  for  multiple  element 
broadside  and  colUnear  dipole  antenna  arrays.  The  multiple  element  concept 
is  studied  as  a  technique  for  eliminating  the  need  for  extremely  high  voltages 
in  the  single  element  type  nuclear  EMP  simulators  by  substituting  the  super¬ 
imposed  electromagnetic  fields  radiated  by  several  lower  voltage  sources  to 
obtain  the  '  ^uivalent  EMP  environment.  E:q)erimentally  measured  input 
impedances  ^.re  compared  with  the  results  from  classical  antenna  theory  to 
determine  the  effects  of  mutual  coupling  between  elements  in  the  array 
coxifigurations.  The  aperture  size  of  the  array  is  shown  to  be  an  important 
factor  in  determining  the  pulse  shape  characteristics  of  the  EMP  environment 
in  the  near  and  far  fields.  Results  of  the  studies  define  specific  advantages 
and  limitations  in  (he  dipole  antenna  array  techniques  for  nuclear  EMP 
simulation. 


t Formerly  with  USA  Mobility  Equipment  R&D  Center,  Ft.  Belvoir,  Vo.  22060 
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1.  INTRODUCTION 


Tho  RES  (Radiating  Electromagnetic  Simulator)  type  antennas  haVe  been  used  success*’ 
fully  as  nuclear  electromagnetic  pulse  (EMP^  simulators.  pThose  simulators  are  dipole  an¬ 
tennas  '  that  are  driven  In  the  time  domain  by  a  MARX  generator  located  at  a  blconlcal  center 
food  section.  A  radiating  dipole  antenna  Is  formed  by  extending  cylindrical  arms  fronl  the 
biconos.  Tho  major  limitation  of  those  types  of  simulators  has  been  tho  rola|tlvely  loW  volt^ 
ages  achieved  to  date  for  a  single  source  pulser.  The  voltages  have  generally  been  less  than 
2  MV  for  operational  dipole  systems.  A  source  in  this  range  cannot  meet  the  large  area 
coverage  requirements  such  as  those  that  exist  for  the  SAFEGUARD  site.  It  is  therefore 
necessary  ,either  I  to  increase  the  voltage  for  a  single  source  or  to  superimpose  the  electro¬ 
magnetic  environments  of  several  sources  to  provide  the  required  threat-level  coverage. 

»  * 

Much  work  has  been  done  In  antenna  theory  on  continuous  wave,  (CW)  dipole  array  sys¬ 
tems.  Such  arrays  have  proved  their  usefulness  In  many  applications.  It  is  of  interest,  there¬ 
fore,  to  examine  the  properties  of  these  arrays  for  application  to  time-domain  EMP  simula¬ 
tion.  In  theory,  one  can  go  from  a  frequency  domain  to  time  domain  by  the  Laplace  or  Fourier 
transformations.  These  are  linear  operations  through  which  some  of  the  useful  properties  of 
CW  arrays  can  be  expected  to  carry  over  to  tlie  time  domain.  . 

The  work  presented  In  this  report  was  conducted  at  the  Electromagnetic  Effects  Labora¬ 
tory  HDL.  Whenever  possible,  theoretical  results  have  been  verified  by  experiment.  Several 
time  domain  dipole  arrays  were  studied  In  order  to  determine  the  possibility  bf  application  of 
such  arrays  to  high  field  level  nuclear  EMP  simulation. 


I 

2.  THE  SINGLE-DIPOLE  ANTENNA  i 

I 

Consider  first  the  single  horizontal  dipole  antenna  to  establish  a  correlation  between 
theoretical  predictions  of  the  radiated  EMP  environments  and  the  experimental  results  ob¬ 
tained  from  tho  models  studied  in  this  report.  The  early  time  EMP  fields  of  a  sin^e-dipole 
antenna  can  be  predicted  with  the  Hertzian  dipole  equations  which,  in  the  frequency  domain, 
are  given  below.  * 


(a.)  =  A  (OS)  c-l'< 


a) 


■>kt|  ^  *0  2 


E,  (oS)  =  h(u)  - 1+ - 1 _ )  cosP 


(2) 


Eg  (OJ)  =  A  (a)  c 


•Jkf  fij±P  +  !£.+ 


sin  0  I 


(3) 


j  c.'  c  r* 


whore  A(oj)  »  I(ai)h/^,  l(aj)  la  the  antenna  current,  h  °  equivalent  length  of  the  antrana, 

Zg  -  /4/«'t  k  -  w/fx  € ,  OJ  -  rod/seo,  /i  and «  are  tho  permoublllty  and  permittivity  of  free  space, 


'Romo,  S.,  Whlnntry,  J.  R.,  Floldi  ond  Wov«  In  Modow  Radio.  Soeend  Edition,  John  WiUy  &'Sont,  Inc.,  N.Y. 
1964,  p.  498. 
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respectively.  The  geometry  for  these  and  subsequently  discussed  com^kinents  is  shown  in 
figure  1,  The  time  domain  expression  for  the  antenna  current  I(t)  is  represented  by: 

I(t)  =  K  \(1  +  a,  t)  e  “*  *  -  (1  +  ttj  t)  e  (4) 


where  K  =  h/47:  determines  the  amplitude  of  the  current,  and  a.^  and  ,  respectively,  de¬ 
termine  the  faU  and  rise  times  of  the  current.  A  plot  of  I(t)  is  shown  in  figiure  2.  The  value 


Figure  1.  Geometry  (or  o  typical  ontenno  element  and  it»  image. 
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Figure  2.  Normalized  input  current  pulse  shape 
computed  from  equation  (4). 


for  the  factor  K  is  determined  experimentally  from  either  a  direct  current  measurement  on 
the  antenna,  or  from  a  measurement  of  in  the  radiation  field  of  the  antenna. 


2.J  Time-Domain  Expressions  for  the  Electromagnetic  Fields 

A  time-domain  solutiot  for  equations  (1),  (2),  and  (3)  is  obtained  by  the  following  procedure, 
llie  variable  is  replaced  by  the  Laplace  variable  s  =  j  ^  and  the  Laplace  transformation  is 
carried  out.  The  results  are  convolved  with  equation  (4)  to  obtain  the  time-domain  expressions 
for  the  dipole  fields: 


where  K  =  h/47r  is  assumed  independent  of  frequency.  This  assumption  requires  that  the  di¬ 
pole  impedance  remains  constant  over  the  pulse  length  of  interest,  and  limits  the  sclution  to 
the  early  time  regime.  That  is,  a  time  <  (2hj/c)  in  figure  1,  where  the  dipole  antenna  does 
not  receive  ground  reflections  which  affect  impedance  of  the  antenna.  In  all  cases  under 
study,  the  antennas  are  located  at  a  height  h,  above  the  ground  so  that  the  rise  time  is  less 
than  the  ground  interaction  time.  We  therefore  expect  a  good  representation  of  the  peak  fields 
from  equations  (5),  (6),  and  (7). 
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Ground  reflections  in  the  radiation  fleld  <Uter  Ihe  free-spare  waveshape  and  must  be  taken 
into  consideration  when  comparing  theoretical  and  measured  results.  In  the  high-frequency 
(early  time)  range,  the  ground  acts  as  a  good  conductor  and  thus  the  theory  of  images  can  be 
applied  by  considering  the  ground  as  a  perfectly  conducting  plane.  The  reflected  (image)  fields 
in  the  presence  of  the  perfectly  conducting  grout'd  plane  have  the  same  form  as  equations  (5), 

(6),  and  (7)  but  arrive  at  the  point  Q  with  a  time  delay  corresponding  to  the  path  length  differ¬ 
ence  transit  time  p/c  between  the  incident  and  image  fields,  where  p  is  the  path  length  differ¬ 
ence  and  c  is  the  velocity  of  propagation.  Orientation  oi  the  incident  and  image  fields  in  figure 
1  assure  that  the  boundary  conditions  are  satisfied  at  the  conductor  surface  (h^  =  0).  We  also 
note  that  the  Fresnel  reflection  coefflcient  assumes  the  correct  limiting  values  of  and  -i  for 
vertical  (electric  field  parallel  to  plane  of  incidence,  horizontal  (electrical  field  per¬ 

pendicular  to  plane  of  incidence,  E,)  polarizations,  respectively. 

The  resultant  fields  in  the  presence  of  the  ground  plane  are  obtained  by  transforming 
equations  (6),  (7),  and  (8)  from  spherical  to  rectangular  coordinates,  and  by  combining  the 
incident  and  image  components  with  proper  time  delay  and  orientation.  This  procedure  gives 
the  X,  y,  z  components  of  the  EMP  environment,  which  corresponds  to  the  measured  fields. 

The  equations  are  given  below,  where  the  "pr  ime”  symbol  denotes  the  imr-ge  fields,  and  the 
subscript  T  represents  the  total,  or  resultant,  field  intensity. 


Hj.  -  5  -  (H^  cos  a  -  cos  I')  y 

Ext ■  cos  'K)  cos  -  (E,  cos  a  -  E'  cos  V)  sin  P 

sina+E'  sin'l')sin'^  -  (E^  sin  a  *  sin  V)  cos 

^  (E^  -  E^)  sin  ,  (E,  -  E;)  cos  tr  (^1) 


2.2  Numerical  and  Experimental  Results 

Equations  (8)  through  (11)  have  been  programmed  on  a  digital  computer  and  numerical 
results  have  been  obtained  for  comparison  with  the  experimental  results.  The  factor  K  in 
equation  (4)  has  been  determined  experimentally  from  the  radiated  to  be  38.42  amperes 
(A).  Figure  3  shows  the  experimental  configuration  of  the  single-dipole  antenna.  A  30'kV 
power  supply  was  used  as  the  source  to  charge  a  pulser,  which  in  turn  fed  a  30-ft  cable  sec¬ 
tion  of  100-n  impedance  to  the  center  section  of  the  antenna.  Three  orthogonal  components  of 
the  electric  and  three  orthogonal  components  oi  the  magnetic  field  were  mapped  at  several 
ranges  on  and  off  the  centerline  (z  =  0  plane)  of  the  mapping  area.  Representative  results  are 
shown  in  figures  4  and  5,  comparing  the  theoretically  and  experimentally  derived  field  com¬ 
ponents.  The  curve  marked  with  the  z  symbol  in  the  magnetic  field  plots  represents  the  free 
space  incident  field  The  other  plots  correspond  to  the  measured  data  shown  in  the  oscil¬ 
loscope  pictures.  The  peak-measu''ed  field  strengths  are  designated  by  each  oscilloscope 
picture  throughout  this  report. 
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SIDE  VIEW 


FRONT  VIEW 
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2.3  Evaluation  of  Results 


The  components  shown  in  figure  4  correspond  to  the  TEM  components  cf  a  plane  wave. 
These  are  designated  PRINCIPAL  COMPONENTS.  Off  the  centerline,  the  environment  be¬ 
comes  noi^lanar  with  the  addition  of  the  E,^  and  the  Ey^.  components.  These  components 
are  not  relatable  to  a  horizontally  polarized  threat  EMP  and  will  be  designated  NON¬ 
PRINCIPAL  COMPONENTS. 

To  correlate  the  simulator  EMP  environment  with  a  particular  threat  EMP  environment 
of  interest,  it  is  necessary  to  know  the  equivalent  free  field  EMP  environment  radiated  by  the 
simulator.  The  EMP  threat  wave-shapes  and  amplitudes  can  then  be  directly  compared  with 
the  simulator  output.  In  general,  the  field-measurement  sensors  and  the  simulator  itself  are 
operated  in  the  proximity  of  the  ground.  Because  of  this,  the  free-space  environment  cannot 
be  measured  directly.  The  theory  of  plane  wave  reflection  at  a  smooth  ground  surface  (or  a 
perfectly  conducting  plane)  can  be  used  to  obtain  the  free  Reid  from  the  measured  tangential 
magnetic  field  (i.e.,  the  magnetic-field  component  parallel  to  the  ground  surface).  The  tan¬ 
gential  component  of  the  incident  magnetic  field  in  figure  1  Is  the  '!x"  component  of  H^.  Boun¬ 
dary  conditions  require  that  the  tangential  component  of  (the  reflected  field)  be  identical 
in  sign  to  the  tangential  component  of  the  incident  field  so  that  the  tangential  magnetic  field 
will  net  vanish  at  the  reflecting  surface.  At  points  above  the  surface,  the  two  fields  add 
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constructively  (since  both  are  of  the  same  polarity)  with  a  time  shift  equal  to  the  path  length 
difference  transit  time  p/c  between  the  incident  and  reflected  fields.  If  the  point  of  field 
measurement  is  close  to  the  ground  surface,  p/c  is  less  than  the  rise  time  of  the  free  space 
incident  {nilse,  and  the  resultant  field  preserves  the  waveshape  of  the  incident  free  field. 

If  such  is  the  case,  the  incident  and  reflected  tangential  components  add  with  very  little  time 
delay,  and  the  only  noticable,  difference  in  waveshape  between  tbe  resultant  field,  and  the 
ftee  field,  is  a  slight  bump  on  the  rising  portion  pf  This  can  be  verified  in  figure  4 

by  comparing  the  measured  and  computed  waveforms  of  with  K^.  We  also  note  that  the 

magnitudes  of  the  computed  and  measured  values  of  agree  closely,  indicating  thht  the  ds- 
sumption  of  perfect  ground  conductivity  is  justified  in  the  early-time  regime  ft  v  the  cony>o- 
nents  that  interfere  constructively  with  ground  reflections.  The  components  that  interfere 
destructively  with  ^  incident  fields  upon  reflection.are  morp  sensitive  to  the  infinite  ground 
conductlvify  assumption,  as  evidenced  by  the  smaller  conqmted  values  of  and  E  com¬ 
pared  with  the  measured  results.  These  componenzs  vanish  at  the  surface  of  a  perfect  con¬ 
ductor,  but  are  grehter  than  zero  at  the  surface  of  an  imperfpet  conductor,  such  as  the  earth. 
Hence,  the  measured  values  .can  be  expected  to  be  greater  than  the  computed  values  near  the 
surface  of  the  groimd.  The  same  trends  appear  for  measurements  of  Hy^  and  at  a  ^si- 
tion  off  the  centerline  as  shown  in  figure  5.  The  measured  wavesh^>s  of  H  ,  however,  difiers 
significantly  off  the  centerline  from  the  computed  pulse  shppe,  pleating  a  discrepancy  between 
theory  and  experiment.  The  difference  is  due  to  radiation  off  the  unterminated  ends  of  the  an¬ 
tenna  used  in  the  experiment,  which  is  not  included  in  the  theore^cal  consideration.  Tbe  adr 
ditional  radiation  results  from  an  abrupt  current  reversal  at  the  end  pf  each  antenna  arm, 
causing  a  magnetio  field  of  opposite  polarity  to  radiate  which  interacts  with  the  principal  mag¬ 
netic 'field  and  changes  the  measured  waveshap^.  The  nocprincipal  electric-field  components 
agree  well  in  pulse  shape  and  amplitude.  Therefore,  for  positions  off  the  centerline,  the  theory 
is  useful  in  predicting  the  relative  amplitude  of  these  conn>onents  compared  with  the  principal 
fields.  The  equations  presented  in  this  section  have  provided;  a  method  for  predicting  tbe  EMP 
environment  radiated  by  a  pulsed-dlpole  antenna  in  tte  presence  of  the  ground  surface.  The 
conqxitations  have  been  suppotted  by  experimental  results.  Ihe  subsequent  sections  of  this 
report  e;ftend  the  methods  developed  for  the  single-dipole  antenna  to  arrays  of  such  antennas 
to  determine  their  possible  value  as  nuclear  EMP  simulators.  , 


3.  FIVE- ELEMENT  BROADSIDE  ARRAY 

,  The  field  strengths  produced  by  a  single-dipole  antenna  nuclear  EMP  simulator  are  limited 
by  the  magnitude  of  the  current  and  voltege  that  can  be  supplied  to  the  antenna  by  a  single- 
pulse  source.  This  section  considers  a  technique  of  increasing  the  simulated  EMP  environ¬ 
ment  by  superimposing  tbe  fields  of  several  simultaneously  activated  antennas  which  are  dis¬ 
tributed  in  an  array  to  form  a  radiating  aperture.  Tbe  five-element  broadside  dipole  array 
model  was  built  to  .study  this  concept.  Illustrations  of  tbe  array  are  given  in  figures  6  and  7. 
The  antennas  are  driven  in  parallel  from  the  30-kV  source  by  five  30-ft-long  100- n  cable  con¬ 
nections.'  This  configurafion  provides,  each  antenna  ydth  a  current  equivalent  to  the  single- 
dipole  current  discussqd  in  section  2.  As  a  result,  the  array  iqmrture  contains  five  times  the 
amount  of  current  compared  to  that  of  tbe  single  sewrro. 


•  3.1  Superposition  of  Electromagnetic  Fields 

,  The  principle  of  superposition  will  be  applied  to  compute  the  EMP  environment  radiated 
by  the  array.  First,  it  is  pointed  out  that  although  five  times  the  single-seuree  ctirrent  emsts 
in  the  array,  the  radiated  environment  is  not  necessarily  assured  to  be  five  times  as  large. 
Ihis  is  a  consequence  of  two  important  factors. 
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Figure  6-  Five-element  broodside  dipole  array  geometry. 


(1)  The  array  is  a  distributed  sourco,  extending  over  several  pulse  risetimes  in  its  di¬ 
mensions.  The  fields  computed  at  some  distant  observation  point  result  Irom  a  retarded  time 
superposition  of  the  individual  antenna  environments.  Hie  retarded  time  refers  to  the  differ¬ 
ent  arrival  time  for  the  EMP  wavefront  from  each  source  to  the  point  of  observation  due  to 
different  path  lengths  between  each  antenna  and  the  observation  point.  Thus,  the  addition  of 
the  fields  from  each  source  is  spread  out  over  a  time  equal  to  the  maximum  transit  time  dif¬ 
ference  between  two  extreme  antennas  of  the  array  and  the  observation  point.  This  results  in 
a  ileld  magnitude  smaller  than  that  produced  a  set  of  simultaneously  superimposed  pulses. 
Consequently,  fields  from  the  five-element  array  are  expected  to  be  less  than  five  times  that 
of  a  single  source  at  finite  ranges  of  observation.  At  an  infinite  distance  of  observation  from 
the  array,  any  finite  aperture  will  appear  as  a  point  source  and  the  retarded  time  efiect  will 
not  be  oteerved.  However,  for  EMP  applications,  energy  limitations  preclude  operation  in 
the  extreme  fim  field  of  the  array. 
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a.  PARALLEL  ARM 
CONFIGURATION 


b.  CABLE  FEED  FROM 
PULSER 

30KV  THRU  30  FT 
of  R69BU  CABLE 


c.  CABLE  FEED  TO 
BICONE  SECTIONS 
OF  EACH  ANTENNA 


(2)  The  second  impoi'tant  factor  to  consider  is  the  possibility  of  losses  occurring  from 
mutual  coupling  between  antenna  elements.  The  radiating  electric  and  magnetic  fields  of  each 
antenna  eventually  strike  every  other  antenna  in  the  array,  and  some  type  of  interaction  (mu¬ 
tual  coupling)  must  occur.  From  elementary  electromagnetic  theory,  one  expects  the  coupling 
to  be  in  the  form  of  induced  currents  which  are  smaller  than,  and  tend  to  oppose,  the  initial 
source  currents  on  the  anteimas.  Thus,  is  addition  to  the  retarded  time  effect,  mutual  coup¬ 
ling  can  be  expected  to  further  reduce  the  magnitude  of  the  EMP  environment  radiated  by  an 
er^rimental  array.  As  a  starting  point  only  the  effects  of  retarded  time  superposition  will 
be  included  in  the  array  field  equations,  and  the  study  of  mutual  coupling  effects  is  considered 
in  a  later  section  of  this  report. 


3.2  Time-Domain  Expressions  for  the  Electromagnetic  Fields 

Ideally,  each  antenna  radiates  the  electromagnetic  fie  ids  described  by  equations  (5),  (6), 
and  (7).  The  equations  corresponding  to  equations  (8)  through  (11)  for  the  single  antenna  are 
expressed  for  the  array  in  terms  of  the  retarded  time  summation  over  the  number  of  antennas, 
N  =  5,  in  the  array: 


N 

n«l 

N 

Myt  = 

n«l 


n 

“n  -  CO'-  COS  COS  a  -  E;^  cos  sin  (14) 


n«l 


Eyx  “n  +  'PJ  sin  sin  +  E^^  sin  cos  9^}  (15) 

n«l 


N 


e.t  -  e;j  cos  e„} 


(16) 


where  the  "prime"  symbol  denotes  the  reflected  (or  image)  field  in  each  equation.  The  n  sub¬ 
scripted  field  components  (xmtribute  to  the  summations  for  the  resultant  fields  only  after  the 
propagation  time  from  the  nth  source  of  the  point  of  observation  Q  has  elapsed.  Figure  1  shows 
the  angles  and  field  components  for  a  ^ical  element  of  the  arri^.  Equations  (12)  through  (16) 
apply  for  any  number  M  of  mutually  parallel  dipole  antennas  which  are  located  at  arbitrary 
positions  parallel  to  the  ground  surface.  The  nth  field  component,  its  corresponding  retarded 
time,  and  the  angle  v^,  and  9^  are  computed  by  using  the  geometry  illustrated  for  a  typi¬ 
cal  source  and  observation  point  in  figure  1.  The  negations  have  been  programmed  on  a  digital 
computer  to  provide  numerical  results. 
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3.3  Numerical  and  Experimental  Results 


Field  measurements  were  conducted  on  the  experimental  model  for  comparison  with  the 
theory.  A  typical  set  of  results  is  shown  in  figures  8  and  9  for  the  observation  point  at  40  ft 
on  and  off  the  centerline  in  the  test  area.  As  was  the  case  for  the  single-dipole  antenna,  the 
computed  plot  represents  the  free  space  radiated  magnetic  field  of  the  array.  The  T  sub¬ 
scripted  components  correspond  to  the  ground  interacted  fields  (l.e.,  the  incident  free  space 
phis  the  giuund-reflected  environments)  which  represent  the  measured  fields  shown  in  the 
oscilloscope  pictures  in  figure  8.  Again  one  observes  that  the  measured  tangential  magnetic 
field  compares  closely  in  waveshape  with  thetheoretical  and  plots,  and  repre¬ 
sents  a  reliable  means  of  determining  experimentally  tlie  free-space  output  of  the  source. 

In  comparing  magnitudes,  the  theoretical  value  for  is  greater  than  the  measured  value. 
This  indicates  that  the  omission  of  mutual  coupling  from  the  computations  has  affected  the 
results.  The  reflection  at  35  nsec  corresponds  to  the  mismatch  at  the  unterminated  ends  of 
the  antennas. 

Figure  9  shows  the  measured  environments  off  the  centerline  of  the  test  area.  In  this 
region,  the  nonprincipal  components  are  significant  in  magnitude,  just  as  they  were  in  the 
single  antenna  case,  and  the  environment  is  nonplanar. 


3.4  Evaluation  of  Results 


A  study  of  the  computed  electric  free  field  indicates  that  the  peak  free  field  for  the  five- 
element  array  is  approximately  five  times  greater  than  the  free  field  of  the  single  antenna. 
Retarded  time  superposition  has  not  affected  the  magnitudes  appreciably  for  the  particular 
array  dimensions  under  study  in  this  section.  This  becomes  apparent  when  the  path  length 
difference  transit  time  is  computed  between  the  extreme  elements  I  and  V  in  figure  6.  The 
time  difierence  is  1.2  nsec  which  is  less  than  the  4  nsec  10  to  90%  free  field  pulse  rise  time 
of  each  antenna.  Consequently,  the  retarded  time  d^adation  occurs  within  a  time  span  of 
1.2  nsec  and  is  beyond  the  resolution  of  the  measurements.  Figure  10  shows  the  computed 
peak  electric  free  field  m^nitudes  for  the  single-  and  five-element  source  as  a  function  of 
range  along  the  center  line. 

The  extension  of  the  sinj^e-dipole  equations  to  a  multiple  soui’ce  array  by  the  method  of 
superposition  has  proved  to  be  a  useful  technique  for  prediction  of  the  array  environment. 

The  conqiutations  have  been  supported  by  the  measured  data.  The  approach  has  demon¬ 
strated  the  an^litude  degradation  due  to  retarded  time  superposition;  also,  it  has  shown  for 
the  particular  array  dimensions  under  study  in  this  section  that  the  effect  is  beyond  the  sen¬ 
sitivity  of  the  measurements.  The  discrepancies  between  the  computed  and  measured  data  are 
therefore  attributed  mainly  to  the  omission  of  mutual  coupling  losses  from  the  computations. 
This  argument  is  supported  by  the  27%  higher  computed  value  of  11,^  over  the  measured 
value  of  tile  component  in  figure  8. 

The  theoretical  and  experimental  approach  has  also  shown  that  the  tree-field  wave- 
shiqw  radiated  by  the  array  can  be  obtained  by  measuring  the  tangential  component  of  the 
magnetic  field,  This  component  interacts  coiistructlvely  with  the  ground  and  preserves 
the  waveshape  characteristics  of  the  free  field.  The  remaining  principal  components  H  and 
E  interact  destructively  with  the  ground  and  therefore  exhibit  the  sharp-spike  waveform 
due  to  interference  from  the  ground-reflected  wave  within  a  few  nanoseconds  after  the  wave- 
firont  arrival  at  the  observation  point. 
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GROUND  RANGE  (ft) 

Figure  10.  Cofflouted  electric  free  field  for  o  single-  ond  five- 
dipole  configuration. 


The  results  of  the  experimental  and  theoretical  studies  discussed  in  this  section  indicate 
that  the  array  technique  can  indeed  be  used  to  increase  field  strengths  radiated  by  EMP  simu¬ 
lators.  The  basic  limitations  of  this  technique  have  been  identified  to  be:  (1)  districted  source 
retarded  time  effects  and  (2)  mutual  coupling  losses.  These  limitations  warrant  further  study 
and  are  discussed  in  the  remalnbig  sections  of  this  report.  .  , 


•  i  I 

4.  COLLINEAR  DIPOLE  ARRAY 

The  coUinear  dipole  array  configuration  shown  in  figures  11  and  12  was  (  ''nstructed  to 
study  the  EMP  environment  of  an  Hrra>  with  a  large  aperture  size  compared  w  ope  rise  time 
of  a  free-field  pulse.  Two  coUinear  banks  containing  three  antenna^  ir  each  form  the  array. 

Each  antenna  has  a  current  ecpiival-ent  to  the  single-dipole  current  discussed  in  section  2.  The 
z  dimension  is  in  the  principal  electric-field  direction  and  corresponds  toi  more  tlum  ten  times 
the  fi:ee  field  pulse  rise  time  of  a  single  element.  In  this  case,  the  different  arrival  times  of 
the  wavefronts  from  antenna  banks  A  and  B  are  clearly  visible  in  figures  13  and  14,  which 
show  data  obtained  at  observation  points  C  and  F  In  the  test  area.  The  initial  peak  in  the 
waveforms  is  due  to  fields  arriving  from  antenna  bank  A.  The  second  peak  is  fields 

from  antenna  bank  B  which  are  displaced  in  time  by  a  delay  correspon^ng  to  the  path  length  dif¬ 
ference  between  the  observation  point  and  points  A  and  B.  The  path  lengths  are  given  in  figurq  11. 
A  comparison  of  transit  time  differences  agrees  closely  with  the  computed  and  measured  results 
in  figures  13  and  14. 

The  field  pattern  of  the  coUinear  dipole  array  differs  distinctly  from  the  field  patterns  of 
the  single  dipole  and  the  broadside  array  configurations.  The  polarity'  of  the  fields  from  each 
bank  of  dipoles  in  the  coUinear  array  is  such  t^t  the  nonprincipal  field  components  (E_^  ,  E,..^, 

)  tend  to  cancel  and  the  principal  field  components  (H^  ,  H  ,  E^  )  tend  to  add  cons&uc-' 
tively  in  the  rectangular  center  region  extending  down  ra^e  beiwem  A  and  B.  Exact  addition 
and  cancellation  of  the  principal  and  nonprincipal  field  components,  respectively,  occur  only 
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TOP  VIEW 


Figure  11.  Collinear  dipole  array  geometry. 

on  the  cenlerllne  of  the  test  area,  defined  along  the  line  DG  in  figure  11.  At  other  points  in  the 
rectangular  center  region,  the  largest  contributions  to  the  BMP  environment  from  each  anten- 
bank  are  due  to  the  6'  independent  and  given  by  equations  (12)  and  (13),  and  to  the 
dominant  sin  0  terms  ,in  equation  (I'O).  These  terms  have  identical  polarity  for  all  antennas  in 
the  array  and  therefore  add  constructively  in  the  field-equation  summations.  Because  of  this 
behavior,  the  colUnear  array  has  the  desired  characteristic  of  reinforcing  principal  field 
conq>onent8  and  reducing  nonprincipal  field  components  within  the  center  region.  Ideally,  this 
should  increase  the  area  of  coverage  by  the  principal  fields  compared  with  the  coverage  by  the 
single  dipole  and  the  broadside  array.  The  results,  however,  show  that  the  gains  in  this  re¬ 
spect  are  relatively  small  for  the  array  dimensions  imder  study. 
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The  ^eatest  limitation  of  the  colUneor  array  is  the  large  time  difference  of  wavefront 
arrival  from  the  two  antenna  banks,  which  pro(luce  the  two  spikes  in  the  waveforms.  To 
eliminate  this  effect,  the  range  of  observation  must  be  increased  to  a  point  adi^re  transit¬ 
time  differences  are  less  than  one  rise  time  of  the  free-iteld  pulse  from  one  antenna.  Delta 
in  figures  13  and  14  illustrate  the  decree  in  translt-tim^  differences  between  the  40-  and 
the  180-ft  ranges.  Even  at  the  far  range,  the  waveshape  degradation  has  not  been  eliminated. 
Points  F  and  G  would  have  to  be  at  least  700  ft  away  from  the  aiTay  to  absorb  retarded  time 
effects  within  a  single-element  pulse  rise  time  at  positions  between  F  and  G. 


5.  EFFECTS  OF  MUTUAi;.  COUPLING 

.  I 

The  previous  sections  predict  the  dipole-array  environments  under  ideal  lossless  condi¬ 
tions.  The  i^)proach  was  successful  in  providing  early  time  waveshape  characteristics  and 
relative  amplitudes  of  the  ^EMP  fields.  One  would  expect  that  the  lossless  case  provides  an 
u{qper  bound  to  the  performance  of  a  dipole  array.  It  is  of  interest,  however,  to  know  the 
extent  of  losses  caused  by  mutual  coupling  between  array  elements.  Such  information  would 
lend  more  confidence  in  equations  (12)  through  (16)  for  use  in  the  design  of  larger  array 
systems.  The  lossless  case  would  give  an  upper  bound  to  the  environment,  and  the  extent 
of  mutual  coupling  would  provide  an  efficiency  factor  for  the  design. 

I 

S.l  Measurements  of  Antenna  Input  Impedance 

I  ‘ 

The  network  representation  of  mutual  coupling  between  antennas  Will  be  used  in  ^hls  study 
The  approach  is  based  on  the  knowledge  of  self  and  mutual  impedance  associated  with  the  an¬ 
tennas.  The  antenna  input  impedance  Z  is  measured  by  time-domain  refledtometry  (TDR). 
Figure  17  shows  the  measured  input  impedances  for  antenna  elements  I,  ID,  and  V  of  the 
broadside  array.  In  each  case,  the  other  four  antennas  weife  removed  for  the  measiurements. 
The  procedure  used  to  make  the  measurement  was  to  transmit  a  step-function  input  pulse 
through  a  50-  D  reference  cable  to  the  antenna,  and  to  observe  the  relative  magnitude  of  the 
reflection  (due  to  impedance  mismatch)  on  a  Tektronix  7704  osciUoscope.  Standard  TDR  pro¬ 
cedure  is  used  to  compute  first  the  reflection  coefficient  -p  at  the  jtmction  of  the  reference 
cable  and  the  dipole  antenna,  and  then  to  obtain  the  antenna  input  iizq)edance  Z  from: 

•  2  -  R  ohms  (17) 

where  R  is  the  reference  5C-n  impedance  and  the  value  of  p  is  giver,  next  to  each  oscilloscope 
picture  in  figures  17  and  18. 

Figure  17  shows  the  measured  dipole  antenna  impedance  as  a  function  of  antenna  height 
above  the  ground.  These  measurements  were  taken  to  test  the  sensitivity  of  the  TDR  tech¬ 
nique.  The  results  show  that  the  ii^xit  impedance  Z  increases  aa  the  antenna  height  above  the 
ground  is  increased.  This  trend  is  correct  if  the  antenna  and  its  image  in  the  ground  are 
considered  as  a  two-wire  transmission  line  with  characteristic  impedance  Zq  given  by  ^ 


^Romo,  S.,  Whinnffy,  J.  R.,  ^i•ldl_and  Wov#s  in  Modern  Rodio.  Second  Edition,  John  Wiley  &  Son*,  Inc.,  N.Y. 
1964,  p.  364. 
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Figure  17.  Time-domain  rellectametry.  Dipole  antenna  impedance  at  a  function 
of  height  above  the  ground. 


Zg  =  120  cosh"*  ,  ^ 


(18) 


where  d  is  the  diameter  of  the  two  conductors  and  s  is  the  separation  distance.  This  equation 
applies  for  a  line  separated  with  an  air  dielectric  and  gives  two  times  the  impedance  value  for 
the  antenna  above  a  lossy  ground  plane.  For  a  perfectly  conducting  ground  plane,  s  is  equal 
to  twice  the  antenna  hei^t  above  the  ground.  The  characteristic  impedance  Zg  can  then  be 
computed  from  equation  (18),  assuming  the  perfectly  conducting  ground  and  an  antenna  diam¬ 
eter  of  4  in.  The  results  are  compared  with  the  measured  data  in  table  I.  It  is  noted  that  the 
computed  values  of  Zg  are  l.S  to  1.8  times  greater  than  the  measured  values  of  7.  The 
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Fijuf*  18.  Time-domain  raflocfomofry.  The  effect  of  odiocent  element*  on  the 
measured  impedance  of  the  dipole  antenna. 


trends  are  corrected,  however,  as  evidenced  by  the  comparison  of  relative  impedance  changes 
as  a  function  of  change  in  antenna  height  for  both  the  computed  and  measured  values  given  in 
table  lb.  The  correlation  of  results  is  good,  adding  to  confidence  in  the  TDR  measurements. 


5.2  Determination  of  Mutual  Coupling  Losses 


To  study  the  cllects  of  mutual  coupling,  the  center  antenna  III  is  chosen  as  a  typical  ele** 
ment  in  the  broadside  array  and  the  network  representation  is  used  to  obtain  the  current  and 
the  electric  field  radiated  by  antenna  in  in  the  presence  of  adjacent  elements  II  and  IV.  The 
mutual  coupling  losses  introduced  by  the  adjacent  array  elements  can  then  be  determined 
comparison  with  performance  of  antenna  in  when  aU  other  antennas  are  remov3d.  Measiure- 
ments  of  impedance  for  antenna  in  were  also  taken  with  aU  five  elements  present  in  the  array. 
No  difference  in  impedance  was  observed  when  elements  I  and  V  were  removed;  therefore,  only 
the  effect  of  adjacent  elements  on  the  performance  of  an  antenna  in  the  array  is  considered. 

Figure  18  shows  the  results  of  TDR  measurements  for  the  input  impedance  of  antenna  Ill 
in  the  presence  of  shorted  and  open  circuited  adjacent  elements.  In  the  open  circuit  condition, 
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Table  la.  Comparison  of  measured  and  com’iuted  input  impedance  for  a  dipole 
antenna  above  a  ground  plane 


Antenna 

Number 

Antenna  Height 
(Meters) 

Computed  Z^ 
(ohms) 

TDK  Measured  Z 
(ohms) 

1 

3.96 

606 

366 

in 

2.44 

548 

320 

_ Y _ 

0.91 

430 

284 

Table  Ib.  Comparison  between  computed  and  measured  relative  impedance 
change  as  a  function  of  change  in  antenna  height 


Confuted 

Measured 

Z(I)/Z(ni)  =  606/548  =  1.11 
Z(ni)/Z(V)  =  548/430  =  1.27 
Z(I)/Z(V)  =  606/430  =  1.41 

366/320  =  1.14 
320/284  =  1.13 
366/284  =  1.29 

the  input  impedance  is  the  same  as  in  the  previous  case  where  no  other  elements  were  present. 
This  agrees  with  discussions  of  mutual  coupling  in  antenna  texts,  namely  that  q;>en-circuited 
elements  should  not  effect  the  performance  of  active  elements  in  an  array.  In  the  short-cir¬ 
cuit  condition,  the  impedance  changes  from  320  to  254  H  indicating  that  tite  adjacent  elements 
have  a  significant  effect  on  the  performance  of  antenna  in.  The  equivalent  circuit  for  this 
condition  is  shown  in  figure  19,  where  e,  represents  the  voltage  at  the  terminals  of  antenna 
in,  and  Zji  and  Zj.  represent  the  self-impedances  of  antenna  in  and  the  adjacent  elements, 
respectively.  The  mutual  inq>edance  is  represented  by  Z^j  .  The  associated  network  relations 
are  given  by  equations  (19)  and  (20). 


0  =  -  Ij  +  Ij  Zjj.  (20) 


Active 

Antenna 


^JJ‘^12  ^22'^12 


Passive 

Shorted 

Adjacent 

Elements 
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The  mutual  impedance  is  determined  experimentally  by  the  following  procedures:^ 

(1)  The  self-impedance  Zjj  of  antenna  in  is  measured  with  all  other  elements  removed. 

(2)  The  impedance  Z  ^  into  antenna  in  is  measured  with  shorted  adjacent  elements  n 
and  IV. 

The  mutual  impedance  is  then  computed  from: 


^12  "  ‘  ^11 


(21) 


Equation  (21)  assumes  adjacent  shorted  elements  with  Z^^  ~  ^22  identical  antenna  ele¬ 
ments.  One  obtains  Z^^  and  Zj  from  figures  17  and  18,  respedtively,  and  the  mutual  impedance 
Zj2  can  be  computed  from  equation  (21): 


Zjj  =  /320  (320  -  254)  =  146  0 


(22) 


Next,  determine  the  voltage  e  1  across  the  biconic  section  of  antenna  ITI  as  follows.  The  gen¬ 
erator  voltage  V,  is  30  kV  at  the  pulser.  This  in  turn  is  transmitted  through  a  100-0  line  to 
the  input  at  the  biconic  section  of  the  antenna.  At  this  point  there  is  a  mismatch  in  impedance 
with  an  associated  reflection  coefficient  F  : 


P  Z-  100 
"  Z  +  100 


(2a 


where  Z  =  320  0,  the  self-impedance  of  antenna  m.  Compute  e.  by: 


e,  =  (l-r)v^  (24) 

In  this  case,  F  =  0.52  and  e,  =  14.4  kV.  Sufficient  data  are  available  to  compute  the  currents 
in  equations  (19)  and  (20).  The  computation  results  in  I,  =  56.83  amperes  (A),  and  Ij  =  -25.9  A. 
Therefore,  in  the  short-circuited  adjacent  element  condition,  antenna  III  induces  a  n^  current 
of  -25.8  A  in  the  adjacent  elements.  This  current  in  turn  radiates  an  opposing  field  to  the  field 
of  the  active  element  ni  and  causes  a  decrease  in  resultant  field  strength. 

A  formula  for  computing  the  peak  free  space  electric-field  strength  at  r  on  the  center- 
line  of  a  bicone-fed  dipole  antenna  is:^ 


60  V  _  60  I 

TT  ~~ 


^ontrni,  R.  C.  Microwov*  Scanning  Antannai,  Voluma  II  Array  Thoory  and  Practica.  Acadamic  Preis,  N.Y. 
pp.  157-160. 

^Schalkunaff,  S.  A.,  Friii,  H.  T.,  Ant»nna»,  Thaory  and  Proctica.  John  Wilay  and  Soni,  Inc.,  N.Y.  1952,  p. 
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(25) 

1966, 

106. 


where  v  is  the  terminal  voltage  and  Z  is  the  self-impedance  of  the  antenna.  In  our  case, 

V  =  Bj  =  14.4  kV,  Z  =  320  n ,  and  r  =  12.2m.  The  net  antenna  current  I  is  equal  to  45A  in  the 
absence  of  other  antenna  elements,  lliis  corresponds  to  a  free-field  strength  of  =  222 
V/m  for  antenna  IE.  In  the  lossy  case,  where  the  adjacent  short-circuited  elements  reduce 
the  net  cxirrunt  to  I  =  56.83  -  25.9  =  30.93A  the  corresponding  field  strength  is  Eg  =  152  V/m. 
(In  the  strict  sense,  the  subtraction  should  include  retarded  time  effects.  However,  the  ele¬ 
ment  spacing  is  so  close  that  a  dii'ect  subtraction  is  justiffed.)  Consequently,  mutual  coupling 
has  reduced  the  efficiency  of  antenna  III  to  69%  of  its  original  output. 

The  mutual  coupling  studies  have  shown  that  significant  losses  can  arise  from  parasitic 
adjacent  elements  in  an  array.  Therefore,  design  considerations  should  include  these  effects. 
The  arguments  presented  in  this  section  have  been  based  on  the  impedance  measurements  and 
the  network  representations  of  mutual  coupling.  From  these  considerations,  we  obtained  the 
net  current  and  field  strengths  required  to  obtain  the  decrease  in  efficiency  of  a  radiating  ele¬ 
ment  in  the  presence  of  adjacent  antennas.  In  the  following  section,  the  efficiency  of  the  arrays 
studied  in  this  report  is  discussed  from  the  standpoint  of  measured  fields.  Thus,  one  has  the 
results  of  two  approaches  to  estimating  array  efficiency:  one  approach  dependent  upon  knowl¬ 
edge  of  impedsmees  and  voltages  in  the  array  (and  independent  of  the  radiated  fields),  and  the 
other  approach  dependent  upon  knowledge  of  the  radiated  fields  (and  independent  of  the  imped¬ 
ances  and  voltages  in  the  array). 


6.  EFFICIENCY  FACTOR  OF  THE  ARRAY  MODELS 


Two  mechanisms  for  degradation  of  performance  have  been  identified  which  introdr 
limitations  in  the  dipole  technique  for  EMP  simulation.  Neglecting  other  losses,  the  physical 
size  of  the  array  determines  a  minimum  range  where  the  waveshape  is  not  degraded  by  the 
arrival  time  differences  from  extreme  elements.  The  array  is  a  distributed  source  for  prac¬ 
tical  ranges:  therefore,  the  field  strengths  do  not  reach  the  ideal  N  times  the  single-element 
output  level,  where  N  is  the  number  of  antennas  in  the  array. 

Define  a  lossless  efficiency  factor  (E.F.)  that  provides  a  measiure  of  bow  much  the  dis¬ 
tributed  nature  of  the  array  affects  the  output,  as  compared  with  the  output  of  a  point  source 
having  N  times  the  single-array  element  strength.  The  H^^  component  is  used  as  a  reference 
field  for  the  conqputations. 


LOSSLESS  E.  F.  = 


'computed  for  lossless  array 
^computed  for  single  element 


X  100% 


(26) 


The  E.F.  factor  approaches  100%  at  ranges  where  the  array  approximates  a  point  source. 
Table  II  shows  computed  values  of  the  lossless  E.F.  for  the  array  configurations  studied  in 
this  report.  The  efficiency  computed  from  is,  in  general,  less  than  that  obtained  from 
the  free-field  output,  as  evidenced  by  the  results  shown  for  the  five-clement  array  in  figure 
10.  For  that  case,  the  ffee-field  efficiency  approached  100%,  whereas  the  use  of  in  equa¬ 
tion  (26)  has  resulted  in  an  efficiency  factor  of  76%. 

Additional  losses  are  introduced  by  mutual  coupling  effects.  These  losses  are  included  in 
a  second  E.F.  which  is  referenced  to  tte  measured  component,  and  which  includes  any 
other  real  losses  present  in  the  array. 
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CONFIGURATION 


Table  n.  The  array  emciency  factor 
REFERENCE  FIELDS  (A/m) 


LOSSLESS 


MEASURED 
H  , 


EFFICIENCY  FACTOR 


LOSSLESS 

MEASURED 

E.F. 

E.F. 

100% 

67% 

100% 

88% 

99% 

76% 

59% 

/i\/  measured  array  \ 

MEASURED  E.  F.  =  ( i )  _ Zl -  x  100% 

\N/\measured  single  element  H  _/ 


Results  of  the  computations  are  given  in  table  11. 

For  antenna  in  with  short-circuited  adjacent  elements,  the  measured  E.F.  is  67%.  This 
agrees  closely  with  the  69%  efficiency  predicted  by  the  mutual  coupling  results.  In  the  open 


I 


circuit  adjacent  element  condition,  the  measured  E.F.  is  88%.  This  indicates  that  losses  are 
introduced  which  were  not  apparent  in  the  impedance  measurements.  The  early  time  fields 
from  the  coUinear  array  can  be  considered  as  a  three-element  broadside  array.  The  meas¬ 
ured  E.F.  is  84%  for  this  condition.  By  comparison,  the  five-element  broadside  array  measured 
efficiency  is  39%.  This  indicates  that  increasing  the  number  of  eijements  in  order  to  increase 
field  strength  has  a  practical  limit,  beyond  which  little  is  gained  by  adding  elements. 


7.  CONCLUSIONS 

I  i 

I 

This  report  summarizes  the  results  of  studies  conducted  at  the  HDL  Electromagnetic  Ef¬ 
fects  Laboratory  on  nuclear  electromagnetic  pulse  simulation  by  dipole  antenna  array  tech¬ 
niques.  It  is  evident  from  these  studies  that  the  aperture  size  and  range  of  operation  are  im¬ 
portant  factors  that  must  be  taken  into  consideration  in  the  dipole  array  design.  The  dimen¬ 
sions  of  interest  are  the  path  length  difference  between  extreme  elements  in  the  array  and  the 
point  of  observation  in  the  test  area.  The  transit-time  differences  must  be  at  le^t  smaller 
than  the  rise  time  of  the  single-element  free  space  pulse  to  prevent  vavesh^pe  degradation. 
Large  apertures  in  either  or  both  the  electric  and  magnetic  directions  may  cause  this  . 
degradation. 

Losses  from  mutual  coupling  must  also  be  considered  In  design  procedures.  The  method 
presented  in  this  report  can  be  used  to  estimate  these  losses  by  taking  impedance  measure-' 
ments  on  reduced  scale  models,  or  possibly  on  full  size  models  of  large  antenna  elements. 

The  efficiencies  computed  for  the  array  configurations  studied  in  this  report  should  be  con¬ 
sidered  as  representative  examples  of  the  correlation  between  the  various  impedance,  current, 
voltage,  and  field  parameters  for  dipole  antenna  arrays.  The  efficiency  factors  indicate  that 
optimal  arrays  should  have  as  few  antenna  elements  as  possible  to  reduce  losses  introduced  by 
mutual  coupling  and  by  distributed  source  effects.  The  theoretical  methods  and  the  computer 
programs  developed  during  this  study  have  been  supported  By  experimental  results,  and  can 
be  used  to  predict  the  output  of  such  arrntfo 
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